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JO  DELING  OF  EQUILIBRIUM  GAS  ADSORPTION  FOR 
MULTICOMPONENT  VAPOR  MIXTURES 

1.  INTRODUCTION 

Single  vapor  equilibrium  adsorption  isotherms  are  often  used  to  assess 

the  relative  efficiencies  with  which  adsorbents,  such  as  activated  carbon, 

remove  specific  vapors  from  air  in  air  purification  schemes.  Single  vapor 

adsorption  isotherms  can  often  be  predicted  at  equilibrium,  subject  to  some 

limitations,  from  the  physical  properties  of  the  adsorbate  vapors  by  tech- 

1  — 

niques  based  upon  the  Dubinin-Pol anyi  concept  of  affinity  coefficient. 

Past  studies  have  provided  vast  amounts  of  adsorption  isotherm  data  for 
single  vapor  adsorption  on  various  adsorbents.  However,  the  practical 
conditions  under  which  adsorbents  are  employed  are  usually  quite  different 
from  the  ideal  laboratory  conditions  under  which  the  single  vapor  isotherms 
are  determined.  For  example,  several  adsorbate  species  are  usually  pre¬ 
sent.  The  vapor  species  to  be  adsorbed  may  not  be  exposed  to  the  adsorbent 
under  conditions  where  equilibrium  can  be  readily  attained.  The  study  of 
mul ticomponent  kinetic  and  equilibrium  adsorption  on  an  adsorbent  is  very 
important;  however,  kinetic  and  equilibrium  adsorption  studies  are  still  at 
a  very  preliminary  stage  while  mixed  vapor  adsorption  studies  are  more 
complex  and  time  consuming.  This  problem  can  be  alleviated  to  some  extent 
if  mixed  adsorption  isotherms  can  be  predicted  from  single  vapor  adsorption 
data. 

Although  there  has  been  considerable  study  involving  the  thermodynamic 
properties  of  adsorbates  on  adsorbents,  relatively  few  studies  have  con¬ 
sidered  adsorption  kinetics.  Studies  to  date^  have  dealt  with  the 
kinetics  of  single  vapor  adsorption  behavior.  Many  of  these  studies  have 
made  use  of  Wheel er’s  approach  to  the  kinetics  of  gas  adsorption  by  beds  of 
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adsorbent  granules.  However,  the  method  has  yet  to  be  applied  to  multicom¬ 
ponent  systems. 

The  development  of  experimental  methods  for  determining  binary  adsorp¬ 
tion  isotherms  and  the  kinetics  of  binary  mixture  adsorption  was  completed 
in  the  first  year  of  the  project.  In  the  past  year,  the  adsorption  of 
CHCl^/CCljj,  n-Hexane/Benzene ,  and  CHgCl2/CHCl^  binary  systems  on  BPL- 
activated  carbon  was  investigated.  BPL  is  a  designation  assigned  by  Calgon 
Corporation.  Results  and  conclusions  from  these  studies  will  be  presented 
in  addition  to  recommendations  for  work  to  be  carried  out  in  future  phase? 
of  the  project. 

2.  BACKGROUND 

2. 1  Equilibrium  Adsorption  of  Mixtures. 

Currently,  thee  are  three  potentially  useful  models  available  which 
have  had  limited  success  in  predicting  the  adsorption  characteristics  of 
mixtures; 

a.  Dubi nin-Polanyi  Pore  Filling  Theory1 

b.  John's  mixture  isotherm  model0 

c.  Myer's  ideal  adsorption  solution  theory^ 


2.1,1  Dubinin-Polany i  Pore  Filling  Theory. 

The  Dubi  nin-Polanyi  theory  has  not  been  used  much  to  predict  multi- 
component  adsorption.  The  results  that  are  available  indicate,  however, 
that  the  theory  has  some  potential  for  application  to  multicomponent 
adsorption.  Bering  and  co-workers1^’11  extended  the  Dubinin-Polanyi  equa¬ 
tion  to  the  adsorption  of  mixtures  by  using  the  following  equation: 


la.  = 

i 


W  9  ZN  log (p  ./p.)  rj 

-  exp  t-BT  (  . . — — )  J 


ZN.v. 

l  l 


(EW 


(1) 
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where  v.  *  the  partial  molar  volume  of  mixture  component  i 

g.  =  the  partial  molar  affinity  coefficient  of  mixture  component  i 
-  the  saturated  vapor  pressure  of  mixture  component  i 
p.  =  the  equilibrium  pressure  of  mixture  component  i 

a.  =  the  number1  of  g  mole  of  component  i  adsorbed  per  gm  of  adsorbent 

»  the  mole  fraction  of  component  i  in  the  adsorbed  phase 

WQ  and  B  =  constants  characterizing  the  adsorbent 

In  practice,  the  quantity  SNjvj  can  be  found  from  the  phase  diagram  of 

the  volume  solution  assuming  a  liquid- like  adsorbate,  while  can  be 

1  2 

found  simply  according  to  an  additive  scheme. 

In  the  theory  of  micropore  filling  in  the  case  of  an  individual 
component,  a  normal  liquid  at  the  given  temperature,  existing  in 
equilibrium  with  its  saturated  vapor  at  the  pressure,  PQ,  is  selected  as 
the  standard  state.  In  the  case  of  multicomponent  adsorption,  it  is 
unclear,  a  priori,  whether  the  state  of  a  solution  whose  composition  is 
equal  to  the  composition  of  the  adsorbed  phase  or  the  state  of  a  solution 
existing  in  equilibrium  with  vapor  whose  composition  is  equal  to  the 
composition  of  the  equilibrium  vapor  above  the  adsorbed  phase  should  be 
selected  as  the  standard  state.  However,  studies  11  have  shown  that 
equation  (1)  is  fulfilled  well  in  both  methods  of  selecting  the  standard 
state.  Selecting  such  standard  states,  we  can  rev/rite  equation  (1)  for  a 
binary  mixture  of  vapors  in  the  following  form: 

W  =  a12V12  =  Wo  “Ptfl  (1°8  (h))2]  (2) 

±2 

where , 


h 


EpSi 


ii 


Equation  2  was  found  to  be  applicable  to  several  systems  J3, 1 4 
2.1.2  John’s  Mixture  Isotherm  Model. 

This  model,  which  was  developed  by  John  and  others,8  assumes  that  the 
single  vapor  isotherms  for  specie?  i  can  be  represented  by  the  following 
equation: 

loglog  P°  =  Ci  +  Di  log  W°  (3) 

where  =  a  constant 
Dj  »  a  constant 

Wj  »  amount  of  adsorbate  in  cm8/g  at  pressure,  p^ 

Pi  -  (P^Psi^  ,qN 

and  superscript  ’o'  denotes  pure  component.  N  =  an  integer  between  2  and 

6. 

A  similar  equation  describes  the  binary  vapor  (components  1  and  2) 
adsorption  isotherm: 

log  log  P1g  =  C12  +  D12  log  W12  (4) 

where  C12  ~  Y-jC.,  +  Y2C2 
D-j  2  =  Y-jD-]  +  Y2D2 

W i 2  “  the  amount  of  mixed  adsorbate 
P1 2  =  (Pl  +  P2)/(Ps1  +  ps2)1qN 

and  and  Y2  are  the  mole  fractions  of  components  1  and  2  in  the  gas 
phase . 

The  constants  C^,  C2,  D-] ,  and  D2  can  be  obtained  from  the  single  vapor 
isotherms  and  then  used  to  calculate  W12,  in  Equation  A,  using  known  or 
assumed  P.  2  values. 

From  the  model,  the  micropore  volume,  WQ,  can  be  computed  for  a  single 
vapor  component  i  as  follows: 

log  WQ  =  (log  log  10"  “  Ci)/Di  (5) 
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Similarly  the  mioropore  volume  in  the  binary  mixture  case  is  given  by: 

log  WQ  =  (log  log  10N  -  C12)/D12  (6) 

John  and  others8  have  shown  that  this  method  can  be  applied  to  binary 

and  ternary  systems  to  compute  their  total  adsorption. 

2.1.3  Ideal  Adsorbed  Solution  Theory  . 

The  ideal  adsorbed  solution  theory  has  been  used  to  predict,  mixed  gas 

q  1  5 

adsorption  with  some  success  at  low  coverages/low  relative  pressure.  ’ 

The  method  assumes  that  the  adsorbed  phase  forms  an  ideal  solution  and 
involves  determining  the  ’spreading  pressure'  for  the  single  vapor  iso¬ 
therms.  The  calculation  is  made  as  follows: 

a.  Obtain  the  single  vapor  isotherms  for  pure  components  in  terms  of 
the  amount  absorbed,  (cc/g)  versus  equilibrium  pressure  (torr), 

b.  The  spreading  pressure  for  these  adsorbates  is  calculated  as 
follows; 


ttA 

RT 


P 

/ 

0 


(7) 


where  it  =  spreading  pressure 

A  =  specific  area  of  adsorbent 

n  =  total  number  of  moles  in  adsorbed  phase/ gm  of  adsorbent 
p°  =  equilibrium  vapor  pressure  of  pure  component 

c.  Calculate  the  vapor  pressure  of  pure  components  at  constant 
spreading  pressure. 

d.  Describe  the  amount  adsorbed  (n^)  at  these  vapor  pressures 
(calculated  in  step  c)  from  step  a. 

e.  Calculate  the  adsorption  equilibria  for  both  components  at  a 
desired  total  pressure,  P,  using  the  following  two  equations: 
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(8) 


py-i  -  p°Ni 

py2  -  p2n2  (9) 


Adding  these  two  equation?,  the  expression  for  an  ideal  liquid  solution  is 
obtained.  p  _  ° 

L  „2_ 

1  "  o  °  (10) 


Pi  -  P2 


The  vapor  phase  composition  is 

Y 


given  by: 


(11) 


N-,  is  the  mole  fraction  of  component  1  in  the  adsorbed  phase  and  is  the 
mole  fraction  of  component  1  in  the  gas  phase. 

f.  The  total  amount  adsorbed  is  obtained  by: 


g.  Finally,  the  amount  of  each  component  adsorbed  from  the  gas 
mixture  is  given  by:  =  nN 

The  complete  isobaric  composition  diagram  for  any  mixture  is  obtained 
by  repeating  the  above  calculation  for  different  values  of  the  spreading 
pressure. 

2.  2  Kinetics  of  Adsorption. 

Studies  on  the  kinetios  of  adsorption  on  activated  carbon  have  been 
previously  reported^-^  for  the  case  of  single  vapor  adsorption  but  have  yet 
to  be  extended  to  multicomponent  mixtures.  However,  many  of  the  mathema¬ 
tical  equations  and  kinetic  processes  which  describe  these  phenomena  for 
single  vapors  should  also  be  applicable  to  mul ti component  mixtures.  Of 
particular  interest  is  the  approach  taken  by  Wheel er^’^  which  has  been 
used  successfully  for  single  vapor  adsorption  kinetics. ^“7  Wheeler's 
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equation,  which  is  baaed  on  the  principle  of  mass  conservation,  can  be 
written  as  follows: 

tb  =  (We/C0Q)[Wb  -  pD  Q  ln(C0/Cx)/kv]  (13) 

where  CQ  =  the  inlet  gas  concentration  in  g/cm-J 

.  i 

kv  =  the  first  order  rate  constant  in  mm 
pB  =  the  bulk  density  of  packing  in  g/cm^ 

W  =  the  kinetic  saturation  capacity  in  g/g 

to 

Wb  =  the  bed  weight  in  gm 
Q  *  the  volume  flow  rate  in  cm^/min 
tb  =  the  breakthrough  time  in  min 
and  Cv  is  the  exit  gas  concentration. 

From  a  plot  of  tb  versus  Wb,  the  saturation  capacity  (We)  and  the 
first  order  rate  constant  (ky)  can  be  obtained.  By  setting  tb  of  Equation 
(13)  equal  to  zero  and  solving  for  Wb  one  obtains 


PR  Q  C 

■  -ir  -  <r>  - 


(14) 


where  is  identified  as  the  critical  bed  weight,  or  that  weight  of  carbon 
just  sufficient  to  reduce  CQ  to  Cx  under  the  test  conditions. 

3.  WORK  OBJECTIVES 


The  ultimate  objectives  of  the  study  are  to  develop  methods  for  pre¬ 
dicting  the  adsorptive  behavior  of  mixed  gas  systems  on  activated  carbon 
adsorbents  from  a  knowledge  of  the  adsorptive  properties  of  the  pure  compo¬ 
nents  and  to  determine  applicability  of  the  Wheeler  equation  to  multi- 
component  kinetic  adsorption. 

The  work  has  been  divided  into  the  following  four  phases:  (1)  deter¬ 
mination  of  equilibrium  adsorption  isotherms  on  BPL-acti vated  carbon  for 
various  binary  vapor  mixtures  at  25^0,  (2)  prediction  of  binary  vapor 
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adsorption  isotherms  on  BPL-actt vated  carbon  for  comparison  with  experi¬ 
mentally  determined  isotherms,  (3)  determination  of  breakthrough  parameters 
through  adsorbent  beds  for  single  vapors  and  binary  vapor  mixtures,  and  (i|) 
testing  the  applicability  of  Wheeler’s  equation  to  binary  systems. 

PROCEDURES  AND  RESULTS 
J4.  1  Binary  Equilibrium  Adsorption 

The  apparatus  and  experimental  procedures  used  for  measuring  the 
single  and  binary  vapor  adsorption  isotherms  have  been  described  in  a 
previous  report  by  Reucroft  and  others.'*® 

The  activated  carbon  was  a  Pi ttsburgh-acfci vated  carbon,  type  BPL,  12- 
30  mesh,  having  an  internal  surface  area  of  about  1000  m2/g  and 
approximately  80$  of  the  internal  surface  area  associated  with  pores  less 
than  2 oft  in  diameter. 

Three  experimental  procedures  were  previously  evaluated  to  determine  a 
precise  method  for  measuring  the  adsorption  characteristic  of  the  binary 
mixtures.  In  methods  A  and  B,  individual  mixtures  containing  various 
relative  pressures  (0.1  0,  0.1  5,  0.20,  0.25,  and  0.30)  of  each  component 
were  made  prior  to  exposure  to  the  carbon. 

In  method  A,  the  same  sample  was  used  for  each  exposure,  with  outgas- 
sing  and  drying  being  employed  betweem  each  exposure.  In  method  B,  a  fresh 
carbon  sample  was  dried  and  outgassed  prior  to  each  exposure.  In  method  C, 
a  mixture  containing  0.3  P/PQ  of  each  component  was  made  and  the  mixture 
was  introduced  as  if  it  were  a  single  vapor.  It  was  concluded  that  Method 
A  was  inaccurate  since  a  residual  amount  of  vapor  could  not  be  outgassed 
under  heating  (300°C)  and  vacuum  (10“5  torr).  In  Method  C,  it  was  assumed 
that  the  gases  mix  ideally  and  behave  like  a  single  vapor  since  the  rela¬ 
tive  volume  of  the  two  gases  adsorbed  is  very  small  compared  to  the  total 
volume  of  gases.  This  assumption  is  valid  since  the  mixed  adsorption 
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isotherms  obtained  by  both  methods  B  and  C  are  in  good  agreement  as  shown 
in  Figure  1.  Since  method  B  is  very  time  consuming  compared  to  method  C, 
we  have  elected  to  use  method  C  for  our  adsorption  studies. 

4.2  Results  and  Discussion 

The  CHCl^,  CCljj  single  vapor,  and  CHCl^/CCl^  mixed  vapor  adsorption 
isotherms  are  shown  in  Figure  2  in  terms  of  amount  adsorbed  (gg-.1 )  versus 
relative  vapor  pressure.  These  mixed  (CHCl^/CCl^)  isotherm  data  were 
obtained  by  Method  C.  As  shown  in  Figure  3,  the  amount  adsorbed  (W-|2,  Sg 
is  a  linear  function  of  chloroform  mole  fraction  in  the  gas  phase.  These 
mixed  isotherms  were  also  predicted  using  the  pore  filling  isotherm  model 


and  John’s  isotherm  model.  The  adsorption  isotherm  data  obtained  from 
experiments  and  the  predicted  isotherm  data  using  the  pore  filling  model 
and  John’s  model  are  listed  in  Tables  1  and  2  and  shown  graphically  in 
Figures  4  and  5.  Although  both  models  underpredict  the  total  amount 
adsorbed  at  higher  pressures,  the  experimental  results  and  the  model  pre¬ 
dictions  are  in  good  agreement  at  lower  pressures. 

It  is  important  to  be  able  to  predict  the  amount  of  each  binary 
component  adsorbed  when  the  amount  of  mixture  adsorbed  (W-j  2)  known.  The 
semi-empirical  function  of  Lewis1  ^  can  be  used,  in  principle,  to  predict 
the  amount  of  each  binary  component  adsorbed: 


W1  W2 


(15) 


where  W-,  +  W2  =  W12  and  W2  and  W12  are  the  amounts  adsorbed  of 

component  1 ,  component  2  and  the  mixture  at  pressures 

o  0  . 

Pi  2  =  p-j  =  P2.  (Note  p-,2  —  P 1  +  P2' 


16) 
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Using  Methods  B  and 


(i-6  6)  M 


Figure  2.  Adsorption  Equilibria  of  Binary  Gas  Mixture  (CHCl^/CCl^ 
BPL-Activated  Carbon  at  25 °C 


on 
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3  -1 

Table  1.  Experimental  and  Predicted  W  (cm  g  )  Values  (CHCl^/CCl^  System) 
Using  John's  Model  and  Pore  Filling  Model  (Mole  Fraction  of 

-3 

CHCl^  =  0.67,  BPL- Activated  Carbon  at  25°C,  Density  =  1.52  g  cm  ) 


3  -1 

Amount  Adsorbed  (W  cm  g  ) 

l  Pi 

Predicted  Values 

Psi 

Experimental 

Join's  Model  Pore 

Filling  Model 

0.0133 

0.3390 

0.3117 

0.3440 

0.0260 

0.3600 

0.3436 

0.3607 

0.034 

0.3680 

0.3530 

0.3668 

0.069 

0.3920 

0.3757 

0.3810 

0.111 

0.4050 

0.3896 

0.3889 

0.149 

0.4140 

0.3977 

0.3931 

0.223 

0.4230 

0.4084 

0.3979 

0.280 

0. 4300 

0.4100 

0.4001 
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3  —1 

Table  2,  Experimental  and  Predicted  W  (cm  g  )  Values  (CHCl^/CCl^  System) 
Using  John’s  Model  and  Pore  Filling  Model  (Mole  Fraction  of 

_3 

CHCl^  =  0-33,  PPL -Activated  Carbon  at  25 °C,  Density  =  1,56  g  cm  ) 


l 

Amount  Adsorbed  (V/  cm  g 

"5 

Experimental 

Predicted 

Values 

P.si 

John's  Model 

Pore 

Filling  Model 

0.002 

0.290 

0.2174 

0.2963 

0,005 

0.315 

0.2826 

0.3229 

0.010 

0.330 

0.3159 

0.3414 

0.020 

0.350 

0.3428 

0.3583 

0.050 

0.378 

0.3725 

0.3772 

0.080 

0.380 

0.3859 

0.3853 

0.100 

0.398 

0.3919 

0.3887 

0.150 

0.412 

0.4023 

0.3942 

0,200 

0.422 

0.4094 

0.3974 

0.250 

0.430* 

0.4147 

0.3996 

^Extrapolated  Value  . 
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CHCI3/CCI4  (Nchc,3=0. 67) 
■  Experimental 
•  Pore  Filling  Model 
▲  John’s  Model 


I  — 

L  Psi 

Figure  4.  Experimental  and  Predicted  Values  of  Total  Adsorption  of 

CHCIq/CCIa  Binary  Vapor  Using  Pore  Filling  and  John  s  Model 
(Mole  Fraction  of  CHCI3  in  Gas  Phase  =  0.67,  Temperature  = 
25°C,  and  BPL- Activated  Carbon) 


Figure  5.  Experimental  and  Predicted  Values  of  Total  Adsorption  of 

CHCI3/CGI4  Binary  Vapor  Using  Pore  Filling  and  John’s  Model 
(Mole  Fraction  of  CHCI3  in  Gas  Phase  =  0.33,  Temperature  = 
25°C,  BPL-Activated  Carbon) 


The  values  of  the  mixture  component?  (CHCl^  and  CC1>,)  calculated 
according  to  the  above  equations  are  marked  by  solid  points  in  Figure  6, 

A?  can  be  seen  from  Figure  6,  the  values  of  CHCl^  and  CCl^  calculated  in 
this  way  show  a  linear  dependence  on  the  mole  fracton  of  CHCl^  In  the  gas 
phase. 

4. 3  Kinetics  of  Adsorption. 

The  experimental  procedure  and  a  schematic  of  the  binary  vapor  test 
apparatus  used  for  determining  the  kinetics  of  adsorption  were  described  in 
an  earlier  report  by  Reucroft  and  others.10  The  concentration  of  each 
component  was  determined  from  the  traced  area  of  the  gas  chromatograph 
peaks.  The  ratio  of  the  exit  concentration,  Cx,  to  the  inlet 
concentration,  C0,  was  plotted  as  a  function  of  time,  t.  The  time  when  the 
breakthrough  concentration  reaches  0.01  CQ  (CX/CQ  =  0.01)  is  defined  as  the 
breakthrough,  t^.  The  exact  breakthrough  time  (tb)  was  determined  as  a 
function  of  carbon  weight  and  composition  of  gas  mixture. 

4.*1  Results  and  Discussion. 

A  typical  breakthrough  curve  for  all  the  single  vapor  runs  shows  a 
sigmoidal  shape  curve.  However,  in  the  case  of  the  binary  vapor  mixture,  a 
vapor  which  is  being  displaced  does  not  follow  this  behavior,  as  evident 
from  the  Cx/C0  ratio  being  greater  than  1.0.  A  linear  regression  analysis 
of  the  tb  versus  Wb  data  (Tables  3-5)  for  several  binary  and  single  vapor 
systems  in  terms  of  the  Wheeler  equation  leads  to  the  determination  of  We, 
W0,  and  ky.  Tables  6-10  show  these  parameters  for  three  binary  systems 
(CHCl^/CCl ||,  CH2Cl2/CHCl2,  and  n-hexane/ benzene).  In  the  case  of  the 
CHCl^/CCl)!  and  CHgClg/CHCl^  binary  systems,  the  experimental  kinetic 
saturation  capacity,  W„ ,  tends  to  show  total  values  which  exceed  those 
obtained  by  the  direct  weight  method  (Wm)  and  equilibrium  gravimetric 
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O  CCI4  (Wheeler’s  Equation) 
•  CCI4  (Lewis's  Equation) 

□  CHCI3  (Wheeler’s  Equation) 
CCI4  *  CHCI3  (Lewis’s  Equation) 
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Binary  Mixture  of  CHCI3/CCI4  (BPL- Activated  Carbon) 


Table  3.  Breakthrough  Time,  as  a  Function  of  Composition  and  Bed 
Weight  (W^)  for  CHCl^/CCl^  Binary  Mixtures* 


Mole  fraction* 
N 

CKC1„ 

-j 

4b 

(min) 

CHC13 

CGI, 

4 

Wb  (g) 

Regression  Equation 

0.00 

— 

6.00 

0.5823 

- 

8.30 

0.7771 

CC.1,  :  t.  =  13.98W  -  2.-18 

— 

9.90 

12.10 

0.8494 

1.0240 

corr.  =  0.9980 

0.20 

5.10 

5.20 

0.  63  77 

CC14:  8.33W  -  0.130 

5.80 

5.95 

0.  7318 

corr.  =  0.999 

6.00 

6.50 

0.7973 

CHC1  :  tb  =  8.681W  -  0.576 

7.30 

7.60 

0.9263 

corr.  =  0.975 

0.35 

4.50 

4.65 

0.5925 

CC14:  9. 903 5 W  -  1.190 

5.90 

6.05 

0, 7466 

fcorr.  =  0.987 

6.85 

7.20 

0.8242 

CHC1  :  9.4210W  -  1.0739 

corr,  =  0,998 

7.45 

7.60 

0.9163 

8.60 

9.00 

1.0250 

0.50 

5.20 

5.70 

0.6108 

CC14:  8.823W  -  0.196 

5.80 

6.00 

0.6810 

corr.  =  0.989 

7.50 

7.90 

0.8511 

CHC1„ :  8.428W  -  0.086 

7.65 

8.20 

0.9191 

corr.  =  0.990 

0.65 

5.90 

6.10 

0.6380 

CC1  :  12.6467W  -  1.8812 

corr.  =  0.949 

6.85 

7.50 

0.7374 

8.80 

10.05 

0.9326 

CHC1  :  10.1801W  -  0.6331 

10.35 

11.60 

1.0743 

corr.  =  0.9997 

0.80 

7.10 

7.60 

0.6856 

CC14:  10.898W  -  0.006 

7.50 

8.60 

0.8111 

corr.  =  0.972 

9.  00 

10.00 

0.8986 

CHC1  :  8.334W  -  1.1930 

corr.  =  0.895 

1.00 

7.90 

_ 

0.6310 

8.80 

- 

0.7520 

CHC1  :  12.05W  -  0.188 

10.  70 

- 

0.8670 

corr.  =  0.9859 

11.70 

— 

0.9868 

*BPL-Activated  Carbon,  Flow  Rate  =  400  cnrVmin,  25°C,  and  P^  =  25  torr. 

total 

^Vapor  Phase 
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Table  4.  Breakthrough  Time,  as  a  Function  of  Composition  and  Bed 
Weight  (Wb)  for  CI^Cl^/CHCl^  Binary  Mixtures* 


Mole  fraction4" 

N 

ch2ci? 

fcb 

(min) 

CH2C12 

chci3 

\  <8> 

Regression 

0.00 

- 

7.00 

10.50 

0. 7000 

1 . 0003 

CHC1_ ;  t, 

3  b 

= 

12.1646W, 

b 

-  1.5738 

- 

12.75 

1.1787 

corr . 

= 

0.9998 

15.50 

1.3994 

0.35 

6.05 

7.50 

0.6998 

C1I  Cl  :  t 

= 

12.1153W, 

-  2.6001 

7.45 

12.00 

10.20 

17-90 

0,8498 

1.1990 

corr. 

= 

0.9975 

CHCV  cb 

ss 

21.0484Wb 

-  7.4179 

corr. 

= 

0.9990 

0.50 

6.00 

8.60 

0.7017 

CH  Ci  :  t 

= 

12.5388W, 

-  3.0496 

4.05 
14 . 70 

15.20 

23.70 

1.000 

1.4007 

Z  Z  D 

corr . 

CHCV  h 

— 

0.9960 

21.5806Wb 

-  6.4839 

corr. 

- 

0.9999 

0.65 

6.00 

10.45 

0.  7004 

CH  C1  :  t 

= 

10.8161W. 

-  1.6632 

8.95 
11 . 90 

16.45 

21.00 

0.9992 

1.2441 

corr. 

CHC1  :  t, 

3  b 

= 

0,9984 

19.4891Wb 

-  2.9901 

corr. 

~ 

0.9971 

1.00 

5.95 

7.55 

- 

1 . 0013 
1.3007 

CH2CV  *b 

= 

7,0211Wb  - 

-  1.2902 

9.55 

- 

1.5505 

corr . 

- 

0.9964 

12.05 

- 

1.8818 

,{BPL- Activated  Carbon,  Flow  "Rate  =  400  cmVmin,  25°C,  and  ^tota2  =  25  torr. 
+Vapor  Phase 
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Table  5.  Breakthrough  Time,  as  a  Function  of  Composition  and  Bed 
Weight  (Wb)  for  n -Hexane/Benzene  Binary  Mixtures* 


Mole  Fraction + 

N  -Hexane 
n 

_ _  ^b 

n-Hexane 

(min) 

Benzene 

"Tb  (s) 

Regression 

0.0 

_ 

5.15 

0.8000 

Benzene : 

- 

6.00 

0.9037 

tk  =  7.7743W.  -  1.0403 

- 

7.55 

1.1006 

- 

9.85 

1.4033 

corr.  -  0.9999 

0.2 

4.50 

4.50 

0.8004 

Benzene:  t^-1 . 447 8Wb~. 0839 

5.95 

5.  95 

1. 0020 

corr.  =  ,9998 

7.40 

7.50 

1.2032 

n-Hexane:  1^=7 . 1996Wt,— 1 . 26: 

corr.  =  1.0000 

0.50 

3.75 

4.10 

0.  7005 

Benzene: 

4.45 

4.50 

0.8515 

t,  =  6.1313W,  -  0,6835 

5.55 

5.85 

1.0018 

6.65 

6.90 

1.2001 

corr.  =  0,998 

8.10 

8.40 

1.4097 

n-Hexane : 

tb  =  6.2459Wb  -  0.5003 

0.80 

3.70 

3.20 

0. 7019 

Benzene: 

4.45 

4.15 

0.8502 

t.  =  10. 11.61  W,  -  4.1374 

7.10 

9-20 

1.3109 

corr.  =  0.9961 

n-Hexane: 

tb  *  5.6232Wb  -  0.2830 

corr.  =  0.9997 

1.00 

4.2 

— 

0.5997 

n-Hexane: 

4.5 

- 

0.6875 

t.  =  8.1474W.  -  1.049 

5.9 

_ 

0.8498 

7.1 

- 

0.9934 

corr.  -  0.9998 

10.35 

- 

1.4028  • 

*BPL-Activated  Carbon,  Flow  Rate  =  400  cmVrain,  25°C,  and  P  =  25  torr. 

total 


Vapor  Phase 
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Table  6.  Equilibrium  and  Dynamic  Adsorption  Data  for  CHCl^/CCl^  Mixture 

on  BPL-Activated  Carbon  (P__  ^  .  =  25  torr.  Flow  Rate  = 

„  n  total 


400  cm\>in  ^ 


and  Temperature  =  25  C) 


Mole  Fraction 

N 

chci3 

(in  Vapor  Phase) 

chci3 

W 

e 

(g/g) 

CC14 

Total  W 

e 

(g/g) 

W 

m 

(g/g) 

W 

s 

(g/g) 

0.00 

- 

0.6962 

0.6962 

0,5674 

0.6281 

0.20 

0.1255 

0.6605 

0. 7860 

0.5560 

0.6240 

0.35 

0.2148 

0.4944 

0. 7092 

0.5799 

0.6200 

0.50 

0.3198 

0.4758 

0.7956 

0.5454 

0.6080 

0.65 

0.4398 

0.3288 

0. 7686 

0.5602 

0.6000 

0.80 

0.5305 

0.1916 

0.7221 

0.5327 

0.5820 

1.00 

0.7815 

- 

0.7813 

0.5047 

0.5700 

Table  7.  Equilibrium  and  Dynamic  Adsorption  Data  for  ri-Hexane/Benzene 
Mixture  on  BPL-Activated  Carbon  (^tota2_  =  25  t0irr  ail<^ 

Flow  Rate  =  400  cm^/min) 


Mole  Fraction 

^n-Hexane 
(in  vapor  phase) 

we 

(g/g) 

Total  Wg 
(g/g) 

W 

m 

(g/g) 

W 

g 

(g/g) 

n-Hexane 

Benzene 

0.0 

- 

0.3300 

0.3300 

0.350 

0.37 

0.2 

0.0666 

0.2501 

0.3167 

0.335 

- 

0.5 

0. 1298 

0.1444 

0.2742 

0.311 

- 

0.8 

0.2040 

0.0894 

0.2974 

0.284 

- 

1.0 

0.2260 

- 

0.2260 

0.250 

0.27 
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Table  8.  Equilibrium  and  Dynamic  Adsorption  Data  for  CH^Cl^/CHCl^ 
Binary  System  on  BPL-Act ivated  Carbon  (P  =25  ton: 

rj  -L  O  t  3  X 

and  Flow  Rate  =  400  cm  /min) 


Mole  Fraction 

N 

CH2C12 

(in  vapor  phase) 

VJe 

(g/s) 

CH2C12  CHC13 

Total  We 
(g/g) 

W 

m 

(g/g) 

0.0 

- 

0.7813 

0.7813 

0.5012 

0.35 

0.193  7 

0.8786 

1.0723 

0.5253 

0.50 

0.2864 

0.6929 

0.9793 

0.5157 

0.65 

0.3212 

0.4380 

0,7592 

0.4993 

1.0 

0.5132 

__ 

0.5132 

0.4517 

Table  9.  Experimental  and  Predicted  Saturation  Capacity  Using  John’s 
and  Pore  Filling  Model  (CHCl^CCl^  Binary  System,  ptotal  = 
25  torr.  Flow  Rate  =  400  cm3 /min,  and  Temperature  =  25°C) 


Mole 

Fraction  * 

N 

CHC13 

W 

g 

(g/g) 

we  (g/g) 

Pore 

Filling 

Model 

John 1 s 
Model 

wm  (g/g) 

wg  (g/g) 

0.00 

0.6962 

0.  6644- 

0.6552 

0.5674 

0.6580 

0.20 

0.7860 

0.6355 

0.6093 

0.5560 

- 

0.35 

0.7092 

0.6284 

0.5997 

0.5799 

0.6130 

0.50 

0.7361 

0.6212 

0.5897 

0.5454 

- 

0.65 

0.7686 

0.6140 

0.5794 

0.5602 

0.6020 

0.80 

0.7221 

0.6068 

0.5687 

0.5327 

- 

1.00 

0.7813 

0.6084 

0.5737 

0.5047 

0.5760 

* 

Vapor  Phase 
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Table  10.  Experimental  and  Predicted  Saturation  Capacity  Values  Using  John’s 
Model  and  Pore  Filling  Model  (n-Hexane/Benzene  Binary  System, 
*Votal  ”  25  torr,  Flow  Rate  =  400  cm3/min,  and  Temperature  -  25°C) 


Wg(g/g) 


Mole  Fraction 

N 

Hexane 

we  (g/g) 

Pore  Filling 
Model 

John  Ts 
Model 

\  (g/g) 

0.0 

0.3300 

- 

0.350 

0.2 

0.3167 

0.3261 

0.3191 

0.335 

0.5 

0.2742 

0.3168 

0.3019 

0.311 

0.8 

0.2994 

0.2918 

0.3260 

0.284 

1.0 

0.226 

_ 

_ 

0.25 
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values  (W„),  However,  in  the  case  of  the  n- hexane/ benzene  binary  system, 

o 

the  W0  values  agree  very  well  with  Wm  and  Wg.  Moreover,  the  We  values  for 
the  n~ hexane/ benzene  and  the  CH2Cl2/CHCl3  systems  generally  Increase  with 
increasing  concentration  of  the  more  volatile  component. 

In  all  three  systems,  it  was  observed  that  the  component  with  higher 
saturated  vapor  pressure  in  a  binary  system  is  displaced  regardless  of  the 
relative  value  of  its  polarity  or  affinity  coefficient.  Also,  the  We 
values  for  all  single  vapors  obtained  from  Wheeler's  equation  at  the  same 
total  pressure  (25  torr)  agree  very  well  with  Wfft  and  Wg  values,  with  the 
exception  of  CHCl^  vapor  (see  Table  11).  It  is  clear  that  the  Wheeler 
equation  is  not  as  effective  when  it  is  applied  to  vapors  such  as  CHCl^  and 
binary  systems  with  CM Cl 3  as  a  component. 

The  breakthrough  time  of  CHCl^  is  faster  than  that  of  CClj,  in 
CHC1  q/CC1  mixtures  but  is  slower  than  that  of  CCl^  in  single  vapor 
determinations.  Similarly,  the  breakthrough  time  of  n-hexane  is  faster 
than  that  of  benzene  in  n-hexane/ benzene  mixtures,  but  is  slower  than  that 
of  benzene  in  single  vapor  determinations.  In  the  breakthrough  studies  on 
CH2C12  and  CHCl^,  CH2C12  displayed  the  faster  breakthrough  time  in  all 
situations . 

The  mole  fractions  of  the  more  volatile  components  in  the  adsorbed 
phase  were  calculated  from  the  We  values.  These  values,  as  a  function  of 
the  mole  fraction  of  the  same  component  in  the  gas  phase  are  shown  in 
Figures  7-9.  From  this  set  of  data,  it  appears  that  the  mole  fraction  of 
CHCl^  in  the  adsorbed  phase  is  similar  in  value  to  the  mole  fraction  of 
CHClg  in  the  gas  phase. 

Figure  7  also  shows  that  the  ideal  adsorbed  solution  theory ^  which 
was  used  previously  to  predict  mixed  gas  adsorption  equilibria  from  single 
vapor  sorption  isotherms  of  the  individual  components,  does  not  predict 
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Table  11.  Kinetic  versus  Equilibrium  Saturated  Capacity  for  Single 

Vapors  at  0.1  Relative  Vapor  Pressure  and  25°C 

3 

(BPL- Activated  Carbon,  Flow  Rate  =  400  cm  /min) 


Saturation 

Capacity  (g/g) 

A 

w 

max 

Vapor 

we  {gig) 

Wg  (g/g) 

(g/g) 

wm  (g/g) 

cci4 

0.696 

0.6281 

0.6678 

0.5674 

CHC13 

0.7813 

0.5700 

0.6586 

0.5047 

ch2ci2 

0.5132 

•  - 

- 

0.4517 

N -Hexane 

0.2260 

0.2680 

0.2805 

0.2500 

Benzene 

0.307 

0.3426 

0.3778 

0.3300 

^Maximum  capacity  of  carbon  (from  equilibrium  gravimetric  studies) 
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Figure  7.  Vapor-Liquid  and  Vapor-Adsorbed  Phase  Equilibria  for  CHCI3/CCI4 
Binary  Systems  (B PL- Activated  Carbon,  400  cm^/min,  25°C) 
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Figure  8.  Vapor-Liquid  and  Vapor-Adsorbed  Phase  Equilibria  for  Benzene/ 
n-Hexane  Binary  Systems  (BPL-Activated  Carbon,  400  cm^/min, 
25°C) 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.< 


Figure  9 


Nch  j  Gig  in  Adsorbed  Phase 


Vapor-Liquid  and  Vapor-Adsorbed  Phase  Equilibria  for 
CHCl^/CHCl^  Binary  Systems  (BPL-Activated  Carbon, 

400  cm3/min,  25°C) 


this.  Thus,  the  ideal  adsorbed  solution  theory  is  not  in  good  agreement 
with  the  experimental  adsorption  data  in  the  CHCl^/CClq  system.  However, 
if  a  liquid  solution  of  CHCl^/CCl^  is  considered,20  the  relationship 
between  the-  vapor  composition  and  the  liquid  solution  composition  is  the 
same  as  that  predicted  from  ideal  adsorbed  solution  theory.  The 
discrepancy  in  the  adsorption  case  may  be  due  to  the  inability  of  the 
Wheeler  equation  to  predict  We  values  in  the  CHClg/CClj,  system.  However, 
in  the  case  of  the  n- hexane/ benzene  system,  normal  behavior  is  observed  from 
20  to  70  mole  percent  of  n-hexane,  i.e.,  the  system  behaves  as  a  normal 
liquid- vapor  equilibrium  (see  Figure  8).  At  lower  than  20  mole  percent  n- 
hexane,  the  curve  becomes  very  flat,  indicating  higher  adsorption  of  n- 
hexane  than  expected  from  normal  liquid- vapor  equilibrium2^.  Figure  8  also 
shows  that  lower  adsorption  of  n-hexane  occurs  at  greater  than  70  mole 
percent  n-hexane. 

In  the  case  of  the  Q-^C^/CHCl^  system,  the  data  covers  a  more  limited 
concentration  range  and  no  conclusions  can  be  drawn  at  the  present  time. 

5.  RECOMMENDATIONS 

Efforts  within  the  last  year  have  deal  t  with  detailed  studies  of 
CHClg/CClj|  binary  and  single  vapors  with  respect  to  equilibrium  adsorption 
isotherms  and  the  kinetics  of  adsorption  on  BPL-acti vated  carbon. 
Additionally,  preliminary  results  have  been  obtained  in  the  cases  of  n~ 
hexane/benzene  and  CH2CI2/CHCI2  binary  and  single  vapors  with  respect  to 
kinetics  of  adsorption. 

For  the  next  phase  of  the  project  several  objectives  have  been 
formulated: 

a.  Continuation  of  binary  vapor  adsorption  studies  on  CHCl^/CCl^,  n- 
hexane/ benzene ,  and  CH2CI2/CHCI2  vapor  mixtures. 
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b.  Extension  of  the  single  vapor  kinetic  adsorption  studies  to  several 
additional  total  pressures  in  order  to  obtain  a  complete  isotherm, 
e.  Further  studies  on  CHCl^  single  vapors  and  mixtures  involving  CHCl^,  in 
order  to  clarify  the  apparent  inapplicability  of  the  Wheeler  equation  to 
vapor  systems  which  contain  CHCl^  as  a  component, 

d.  Application  of  the  pore  filling  model,  John*s  model,  and  the,  ideal 
adsorbed  solution  theory  to  the  n-hexane/benzene  and  CH2Cl2/CHClg  binary 
systems . 
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